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Abstract 

I review the deconfining phase transition in an SU{N) gauge theory 
without quarks. After computing the interface tension between Z{N) 
degenerate vacua deep in the deconfined phase, I foUow Giovannangeh 
and Korthals Altes, and suggest a new model for (discrete) Polyakov 
loop spins. Effective theories for (continuous) Polyakov loop spins are 
constructed, including those with Z{N) charge greater than one, and 
compared with Lattice data. About the deconfining transition, the ex- 
pectation values of Z{N) singlet fields ( "quarkless baryons" ) may change 
markedly. Speculations include: a possible duality between Polyakov loop 
and ordinary spins in four dimensions, and how Z{N) bubbles might be 
guaranteed to have positive pressure. 

To appear in the Proceedings of "QCD Perspectives on Hot and Dense 
Matter", a NATO Advanced Study Institute, August, 2001, Cargese, France 

1 Overview 

In these lectures I review the deconfining phase transition in a "pure" SU{N) 
gauge theory, without dynamical quarks. Gauge theories are ubiquitous in 
physics, so their phase transitions are manifestly of fundamental importance. 
Two examples may include the collisions of large nuclei at high energy and the 
early universe. 

The phase transitions of gauge theories without quarks are of especial inter- 
est, since then the order parameter, and many other aspects of the phase tran- 
sition, can be characterized precisely |[ ||]. While of course QCD includes 
quarks, this is not an academic exercise. Recent results from the Lattice on 
"flavor independence" — for both the pressure ^, ||] and quark susceptibilities 
[|j — suggest that the results from the pure glue theory may be, in a surprising 
and unexpected fashion, relevant for QCD. (Whether flavor independence can 
be generalized when there are many light flavors is not known.) 

Albeit indirectly, the Lattice has already told us much about what happens 
in a pure gauge theory with three colors. By asymptotic freedom, at infinite 
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temperature the pressure is that for an ideal gas of gluons. In the confined phase, 
the pressure is very small, essentially zero. So the question is, as the pressure 
turns on at the transition temperature Tc, how rapidly does it approach the ideal 
gas limit? The Lattice tells us relatively quickly: by about 2Tc, it is already 
80% of the way to ideality. The "2" in 2Tc is meant schematically; it is certainly 
not, say, lOTc. Above 2Tc, the pressure then approaches ideality slowly, from 
below. 

This suggests that from temperatures of 2Tc on up, that the theory is some 
sort of quasiparticle gas of thermal quarks and gluons; i.e., a Quark-Gluon 
Plasma. By this I mean that after suitable dressing from bare into quasi- 
particles, the residual interactions are weak. This is seen from the Lattice: 
like the ratio of the true to the ideal pressure, the ratio of (gauge-invariant) 
masses to the temperature also vary slowly above 2Tc 0, ^ . 

It is known that for the free energy, direct perturbative calculations fail at 
astronomically high temperatures B; there is only a perturbative Quark-Gluon 
Plasma above temperatures of ~ 1(F GeV. Thus from temperatures of 10^ GeV 
down to 2Tc, the theory is what I call a non-perturbative Quark-Gluon Plasma. 

One approach to the non-perturbative QGP is to fold in the effects of Debye 
screening. It is known that the pressure, as obtained from the Lattice, can be 
fit to an ideal gas of massive gluons all the way down to Tc . The problem 
is that as introduced, these masses aren't gauge invariant. Further, in the end 
one is just fitting one function of temperature, the pressure, to another, thermal 
masses. (Still, it is most intriguing that these fit thermal masses become large 
near Tc.) 

Another approach is provided by the resummation of Hard Thermal Loops 
(HTL) |12|, For the gauge field, if Aq is the time-like component of 

the vector potential, the Debye mass term is just ^ tr(A§). For gluons, HTL's 
are the gauge invariant, analytic continuation of this mass term from imaginary 
to real time. At present, though, HTL resummation has been used mainly to 
compute the pressure. The crucial test, yet remaining, are the results which it 
gives for Polyakov loop correlation functions ||^, |l^] . 

A method to compute all static correlation functions in the non-perturbative 
Quark-Gluon Plasma starts with the (perturbative) construction of an effective 
theory in three dimensions |l6| . From the original gauge theory at T ^ 0, 
static magnetic fields produce three dimensional gluons, while static electric 
fields give Aq, as an adjoint scalar coupled to these gluons. Due to the power- 
like infrared divergences of gauge theories in three dimensions, perturbative 
calculations are useless in this effective theory. Since the effective theory is 
purely bosonic, though, static correlation functions can be efficiently computed 
by numerical means on the Lattice. While heroic perturbative calculations are 
required, for three colors this method appears to work down to 2Tc, where the 
method itself indicates its failure. To be fair, this method only yields static 
correlation functions, and not those in real time. With HTL resummation, this 
continuation is almost automatic. 
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What about below 2Tc? If the transition is strongly first order, then presum- 
ably the quasiparticle regime extends all of the way down to Tc- Lattice data 
suggests that for four colors the transition is strongly first order. (Although 
data on correlation lengths near Tc is absent.) For more than four colors, if the 
deconfining transition is like the Potts model, then it becomes more strongly 
first order as N increases. On the Lattice, "reduced" models have a first order 
transition at infinite N [p^ . In the continuum, for really no good reason, I 
remain unconvinced [p^ . 

For two colors, however, the deconfining transition is almost certainly of sec- 
ond order |2^. For three colors, the transition is of first order [|T), as predicted 
by Svetitsky and Yaffe Recent results suggest, however, that the transition 
is so weakly first order that it is more accurate to speak of a "nearly second 
order" transition |Q. 

This is seen most clearly not from the pressure, but from the behavior of 
electric and magnetic masses. I define the electric, niei, and magnetic masses, 
Tnmag, in ^ gaugc invariant way, from the fall-off of the two-point functions for 
Polyakov loops, and the (trace of the) magnetic field squared, respectively. In 
the perturbative QGP, mei > rumag 0, while in the non-perturbative QGP, 
rUei ~ rrimag |§| ■ In contrast, in the transition regime rrimag/T is approxi- 
mately constant, but rUei/T appears to drop by a factor of ten as T : 2Tc r+ 
[ p2| . A similar drop in the string tension is seen as T : ^ . 

The broad outlines of the appropriate effective theory in the transition region 
for a (nearly) second order deconfining transition have been known for some time 
^, |[ |2^, In the QGP regime, one deals with Aq. In the transition region, 
one trades Aq in for the thermal Wilson line. Traces of powers of the thermal 
Wilson line give Polyakov loops, which are gauge invariant. 

The effective theory of Polyakov loops is just beginning ||2^, ^ |2^, 

0. If true, instead of the Aq quasiparticles of the QGP, near Tc it is more 
useful to view the theory as a condensate of Polyakov loops. This could produce 
dramatic signatures in heavy ion collisions, with hadronization at Tc computed 
semiclassically from the decay of Polyakov loop condensates H^, |2^, . 

In these Lectures I provide some background to understand these questions. 
After explaining why SU{N) gauge theories have Z{N) degenerate vacua B, 

1 compute the interface tension between these vacua at high temperature [ p2| , 
^ lU, 11^, ^ m, Viewing the Z{N) vacua as discrete spins, and 
using results of Giovannangeli and Korthals-Altes I propose a new spin 
model, distinct from the usual Potts model. (The order of the transition seems 
to agree with Potts for A^ < 4, but is unknown for A^ > 4.) I then develop 
an effective theory of Polyakov loop spins, considered as continuous variables. 
Even with the limited amount of relevant Lattice data which exists at present, 
the form of this effective theory is sharply constrained. Especially intriguing is 
the possibility that the expectation values of fields which are singlets under the 
Z{N) symmetry — which I term quarkless baryons — change suddenly about 
Tc. 
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The lectures are in part pedagogical, in part base speculation. The latter 
includes a possible duality in four dimensions, and how Z{N) bubbles, which in 
perturbation theory can have negative pressure [ p9[ , might be ensured of positive 
pressure non-perturbatively. I also review what is known about renormalization 
of Wilson loops and Wilson lines ^ ||, ||. 

In these lectures I only discuss the Polyakov Loops Model ^ |2^, 

|3l| in passing. I hope to provide a basis for understanding its motivation. At 
present, a major unsolved problem is the analytic continuation from imaginary 
to real time. In abelian gauge theories, the analytic continuation of the Debye 
mass term gives the Random Phase Approximation p^ ; in nonabelian theories, 
it gives Hard Thermal Loops Absent any results whatsoever, I do have the 
temerity to coin a phrase for the analytic continuation of the Polyakov Loops 
Model to real time, as the Nonabelian Random Phase Approximation^^ . Time 
(dependence) will tell. 



1.1 Saturation 

My motivation for studying this subject is its relevance for the collisions of large 
nuclei at very high energies. Thus at the outset, I wish to add some general 
comments about nucleus-nucleus collisions, and especially how it might relate 
to models of saturation. 

Consider a completely implausible situation, the collision of two neutron 
stars (say) at relativistic energies. For a neutron star, the transverse area is es- 
sentially infinite on nuclear scales. If the stars completely overlap (zero impact 
parameter), then at very high energies, a nearly baryon-free region is generated 
between them (about zero rapidity). The system starts out with energy den- 
sity, but no pressure. It is then reasonable to think that the system builds up 
pressure, and thermalizes, before it flies apart. 

The crucial question for the collisions of two large nuclei is whether for gold 
or lead nuclei, with A ~ 200, that this finite value of A is close to infinity, 
or represents some intermediate regime. This will be decided by experiment, 
at the SPS, RHIC, and the LHC. After one year of running, it appears clear 
that something dramatic has happened between SPS and RHIC energies [ ^Sf . 
Precisely what is still being sorted out. 

I wish to comment here on the relevance of "saturation" at these energies 

. At a very pedestrian level, this can be viewed as a type of finite size effect 
at A < cxD. Consider a collision in the rest frame of one nuclei. The diameter 
of the other nucleus, with A ^ 200, is ~ 15 fm. This distance becomes Lorentz 
contracted. Thus we can ask, in the rest frame of one nuclei, when does the 
incident nuclei look like a pancake of negligible width? We want the distance 
to be really small on typical hadronic scales. If a typical hadronic scale is 1 fm, 
then, a small scale might be 1/4 ^ 1/3 fm, say. To contract 15 fm down to these 
sizes then requires a center of mass energy per nucleon pair, ^/s/A, on the order 
of ^ (3 ^ 4) X 15 — 45 — > 60 GeV. While an extremely naive estimate, this 
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does seem to be a reasonable estimate for the \/s/A where the distribution of 
particles in AA collisions changes dramatically, developing a "central plateau" 
as a function of pseudo-rapidity. Thus perhaps the appearance of the central 
plateau is where the effects of saturation first appear. 

The details are far more involved, but but our understanding of saturation, 
which is known formally as the Color Glass [Q, gives us some confidence that 
the system is described in terms of "saturated" gluons, with a characteristic 
momentum p^at ^1 — 2 GeV. 

The Color Glass changes all assumptions in one fundamental respect. Fol- 
lowing Bjorken, the usual assumption is that the system thermalizes on a "typ- 
ical" hadronic time scale, ~ 1 fm/c. If saturation kicks in, however, all typical 
scales are then given in terms of 1/psat, which is a much smaller time scale, 
~ .1 — .2 fm/c. If this is the relevant scale, then evolution to a thermal state, for 
a nucleus of size 6 fm, appears much more plausible. Certainly it yields testable 
predictions, which are testable experimentally. 



2 Z{N) symmetries in SU{N) 

I start by reviewing how, following 't Hooft a global Z{N) symmetry emerges 
from a local SU (N) gauge theory. The action, including quarks, is 

C=^trGl, + qipq, (1) 

where 

D^. = df,- igAf, , = — [Df„D,] ; (2) 

= A^^t", with the generators of SU{N) normahzed as tr(t"t'') = (5"V2. The 
Lagrangian is invariant under SU{N) gauge transformations Q, 

D^-^n^D^n , q-^n^. (3) 

As an element of SU{N), fl satisfies 

n^n = 1 , detf7 = 1 . (4) 

Here £7, as a local gauge transformation, is a function of space-time. 

There is one especially simple gauge transformation — a constant phase 
times the unit matrix: 

= e"^l. (5) 
To be an element of SU{N), the determinant must be one, which requires 

= ^ , j = 0,l...(7V-l). (6) 

Since an integer cannot change continuously from point to point, this defines a 
global Z{N) symmetry. 
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2.1 Z{N) at nonzero temperature 

As a particular gauge transformation, Z{N) rotations are always a symmetry 
of the Lagrangian, either with or without quarks. I now show that with quarks, 
they are not a symmetry of the theory, because they violate the requisite bound- 
ary conditions. 

I work in Euclidean space-time at a temperature T, so the imaginary time 
coordinate r, is of finite extent, r : — > /? = 1/T. The proper boundary 
conditions in imaginary time are dictated by the quantum statistics which the 
fields must satisfy. As bosons, gluons must be periodic in r; as fermions, quarks 
must be anti-periodic: 

A^ix,P) = +A^ix,0) , qix,(3) = -q{x,0) . (7) 

Obviously any gauge transformation which is periodic in r respects these bound- 
ary conditions, 't Hooft noticed, however, that one can consider more general 
gauge transformations which are only periodic up to Qc- 

n{x,(3) = n^ , n{x,o) = i. (8) 

Color adjoint fields are invariant under this transformation, while those in the 
fundamental representation arc not: 

= QlAf,{x,P)n, = A^{x,l3) = +A^{x,0), (9) 

g"(f,/3) = nU{x,l3) = e-^^(x,/3) ^ - q{x,0) . (10) 

Here I have used the fact that flc, as a constant phase times the unit matrix, 
commutes with any SU{N) matrix. Consequently, pure SU{N) gauge theories 
have a global Z{N) symmetry which is spoiled by the addition of dynamical 
quarks. 

In the pure glue theory, an order parameter for the Z{N) symmetry is con- 
structed using the thermal Wilson line: 

L(x) = P exp j Ao(x, r)dr^ ; (11) 

g is the gauge coupling constant, and Aq the vector potential in the time di- 
rection. The symbol P denotes path ordering, so that the thermal Wilson line 
transforms like an adjoint field under local SU (N) gauge transformations: 

L(f) ^ n^^ix, f3) L{x) n\x, 0) . (12) 

The Polyakov loop is the trace of the thermal Wilson line, and is then gauge 
invariant: 

^=ltrL. (13) 
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Under a global Z{N) transformation, the Polyakov loop £i transforms as a field 
with charge one: 

£ e"''e . (14) 

At very high temperature, the theory is nearly ideal, so g « 0, and naively 
one expects that {£) 1. Instead, the allowed vacua exhibit a A^-fold degener- 
acy: 

{£) = exp(^^^ eo , J ^ 0,l...(iV-l), (15) 

defining to be real; £o 1 as T — > oo. Any value of j is equally good, and 
signals the spontaneous breakdown of the global Z{N) symmetry. 

At zero temperature, confinement implies that £q vanishes The modulus, 
£o, is nonzero above Tc'. 

£a = , T<T, ; £o > , T > . (16) 

As is standard, if £o turns on continuously at Tc, the transition is of second order; 
if it jumps at Tc, it is of first order. What is atypical is that the Z{N) symmetry 
is broken at high, instead of low, temperatures. For a heuristic explanation of 
this in terms of Z{N) spins, see sec. (3.2). 

One can also understand what it means to say that the global Z{N) sym- 
metry is violated by the presence of dynamical quarks. In the pure glue theory, 
Z{N) rotations take us from one degenerate vacua to another, all of which have 
the same pressure; see sec. (3.2). Adding dynamical quarks (with real masses), 
the stable vacuum is that for which {£) is real, j — Q As discussed 

in sec. (3.4), because of quarks the pressure for a Z{N) state with j 7^ is 
less than the stable vacuum; thermodynamically, they are unstable. What is 
exciting to some of us [Q, however, is the possibility that these Z{N) rotated 
states might be metastable. Such "Z{N) bubbles" could have cosmologically 
interesting consequences [Q. Others contest whether any of this makes any 
sense 

The usual interpretation of the Polyakov loop is as the free energy of an 
infinitely heavy test quark |47| : 

{£) =e^p{-Fte,t/T) . (17) 

This cannot be quite right: when N — 2, the left hand side can be of either 
sign, while for > 3, it is complex. In contrast, free energies are real, so the 
right hand side is positive. 

I suggest a different view. Consider the propagator for a scalar in a back- 
ground gauge field (the extension to fermions is irrelevant here). This propaga- 
tor is given by a Feynman sum over paths: 

Vx'^exp f- [ ds(— + m + igA'^xA^ , (18) 
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with s the length of the path for the worldhne of the particle, and x = dx^/ds. 
A very heavy quark moves in a straight line; in imaginary time, it sits wherever 
you put it. As a colored field, however, it also carries a color Aharonov-Bohm 
phase. This phase is nontrivial, and is precisely the thermal Wilson line. Thus: 
the Polyakov loop, £, is the trace of the propagator for a test quark. 

Confinement then means that (the trace of) this propagator vanishes. For 
two colors, for example, the confining vacuum is Z{2) symmetric: it is composed 
of domains, of definite size, in which i = +1 and £ = —1. As the test quark 
travels through each domain, it picks up one phase or the other. Over a very 
long path, these phases cancel out, giving zero overall. The same holds for 
higher A^, except that there are then N types of domains. This picture appears 
analogous to the localization of an electron in a random potential |Q . 

Z{N) rotations can be expressed in terms of the canonical formalism H, |35| . 
In = gauge, usually the partition function is strictly a trace of exp{—H/T) 
(TL is the corresponding Hamiltonian) , sandwiched between the same state: 

Z = I](V^|exp(-W/T)|^). (19) 

The sum is over all gauge invariant states in the theory. In a canonical for- 
malism, one inserts a projector to ensure that the states are gauge invariant, 
satisfying Gauss' Law, although the process is standard. Instead, what enters 
here is a "twisted" trace: 

Z(f},) = S](V''lexp(-W/r)|V'). (20) 

Here, represents the gauge transform of -0 by the gauge transformation fie- 
This is not automatically equal to -0, because Gauss' Law only ensures that 
the state is invariant under local gauge transformations, and does not restrict 
its behavior under the global gauge transformations. This twisted trace is only 
possible in a gauge theory. 



3 Z{N) Vacua and Bubbles 

3.1 Tunneling between Degenerate Z{N) Vacua 

Typically, discussions of gauge theories at nonzero temperature work up from 
zero temperature. But the zero temperature theory confines, which is compli- 
cated. Instead, I work down from infinite temperature, in a gas of nearly ideal 
gluons. I now compute the amplitude to tunnel from one Z{N) vacua to another 
by semi-classical means U, ^ . 

For simplicity, I work with two colors. To compute the interface tension, I 
put the system in a box: 

t:0^/3 , x,y:0-^ Lt , z : ^ L . (21) 
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I choose one arbitary direction, say the z-direction, and make that much longer 
than the transverse spatial directions, x and y, and than the direction in imag- 
inary time, T. I impose boundary conditions such that the Polyakov loop has 
one value in Z(2) at one end of the box, and the other value at the other end 
of the box: 

£(0) = 1 , £{L)=-l. (22) 

With these boundary conditions, the system is forced to form an interface be- 
tween the two ends of the box. The simplest interface is one which is constant in 
the transverse directions. Thus the natural expectation is that with the above 
boundary conditions, the action is 



Sinter = , C = Cq + QCi + . . . (23) 



One expects the result to start as ~ as a semiclassical probability for 

tunneling in weak coupling. This is standard with instantons, etc. Higher order 
corrections to the leading term, co, are generated by including quantum effects, 
and produce the corrections ci, etc. 

As an ansatze for the Z{N) interface, I take 

<(z) = ^q{z)as , ^3 = ( J ) (24) 

With this ansatz, the Polyakov loop is 

i{z) = cos {nq{z)) . (25) 
Thus the boundary conditions are satisfied by taking 

(7(0) = , g(L) = 1 . (26) 
At the classical level, the action of the above configuration is: 

r0 



= dz . (27) 



Unsurprisingly, the action for the gauge field becomes a kinetic term for the 
classical field. There is only a kinetic term, since the classical field commutes 
with itself. 

This implies, however, that at the classical level, there is no difference be- 
tween the two vacua, or indeed any state with q 0\ One can take q(z) = z/L; 
then the action is ~ L/L'^ and vanishes as L ^ oo. In terms of the 

above, 

Co = . (28) 
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I now show that quantum corrections generate a nonzero value for ci . Since 
Co vanishes, this is then the leading term in a semiclassical expansion; the tun- 
neling probability is then not ^ l/g^, but only ^ 1/g. 

To show this, it is necessary to compute the quantum corrections about the 
above semiclassical configuration, taking 

A, = + Al . (29) 

The computation is a bit involved, but an excellent exercise in the use of the 
background field method |^ , which is always good to know. 
It is convenient to take background field gauge, 

D^jAl = d,Al - = 0. (30) 

I work in euclidean space-time with (+ + H — h) metric. Note that the appropriate 
covariant derivative is that in the adjoint representation. With this gauge fixing, 
the Lagriangian density is 

£ = i tr (G^f + 1 tr {D^^ ■ A'^f + r] {-D^' ■ D) , (31) 

suppressing ugly vector indices. 

With this form, it is easy integrating out the quantum fields to one loop 
order, and obtain the quantum action 

S'{A'') - itrlog(A;J) - trlog(A-i) . (32) 
At one loop order, the full effective action is the sum of the classical action. 



(P7|), and the quantum action, (32). The quantum action involves the inverse 



propagators in a background field. That for the gluon is 

A;.^ = -Dl5^, + (1 - r')D''^D:' + 2zg [G^'„ ] , (33) 
while that for the ghost is (to lowest order in g) 

= -Dl, . (34) 

These results are valid for an arbitrary background gauge field. 

I now make a crucial assumption, and assume that the field q{z) is constant 
in space. For the relevant tunneling amplitude, in fact q{z) does depend upon 
z; what happens, however, is that for the quantum action, this variation only 
enters to higher order in the coupling constant. 

This assumption vastly simplifies the problem. Since Aq lies in the 173 direc- 
tion, it is a diagonal matrix, and covariant derivatives commute: 

[D^^D:'] - G^^. = 0- (35) 
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For example, one can easily show that the quantum action is independent of 
the gauge fixing condition. The variation of the quantum action with respect 
to the gauge fixing parameter ^ is 



d 



I tr {-D'^D'ijA^) , (36) 



with A^'j, the gluon propagator. Normally, this is difficult to compute. In the 
present example, however, if covariant derivatives can be assumed to commute 
with each other, then they can be treated just like ordinary derivates, so that 

Consequently, 

^^" = 1^1). (38) 

Thus there is gauge dependence in the quantum action, but it is completely 
independent of the background field, and so can be safely ignored. 
Consequently, I take background Feynman gauge, ^ = 1, and 

S<^ = tr log{-Dli) . (39) 

The overall factor is expected for a massless gauge field, with two (spin) degrees 
of freedom. 

To compute the determinant in this background field, I introduce the "lad- 
der" basis. 

This is useful because of the commutation relations: 

[(73, a±] = ±2<7± , (41) 
so that the covariant derivative becomes 

D^'a^ = (^do - ig (^^q^ [(73, ]^a^ = I (2t,T) (n ± q)a^ . (42) 

In the last expression, I have gone to momentum space. Remember that given 
the periodic boundary conditions at nonzero temperature, 

ko = i 2'kT n , (43) 

where for a bosonic field, such as a gluon, the periodic boundary conditions 

require that n be an integer, n = 0, ±1, ±2 In the present case, it is handy 

to introduce the shifted momentum, 

= i 2'kT {n±q) . (44) 
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In the trace, the sign of q doesn't matter, so that in momentum space, 

S" = 2 tr log + . (45) 

In computing integrals at nonzero temperature, the usual approach is to do 

the sum over the fco's first by contour integration or the like, and then integrate 
over the spatial momentum. In the present example, instead it is better to first 
integrate over the spatial momentum, and then sum over the fco's, using zeta 
functions tricks. 

Only the variation of the action with respect to q, 

a _ o..^t^ f cPk h 



is needed. The integral over k can be done using dimensional regularization, 
viewing like a mass. Using the standard integral. 



/ 



d-k ^ = na-n/2) TT"/^ 

(fc2+rn2)« r(o) (to2)»-»/2 ' ^ ' 



the result is finite: 

fi / 1 

n= — cso ^ 

+ OC 

= -S-k'^T^LIL y {n + q)\n + q\. (48) 

n= — OO 

The sum over n, where n runs from minus infinity to plus infinity, is turned into 
a sum from zero to plus infinity: 

= -Stt'T^^L ^ ((n + qf - (n + (1 - q)f) (49) 

While these sums are very divergent, mathematicians know how to handle them. 
They are defined by the analytic continuation of the Riemann zeta-function: 

Using 

C(-2,«) = ' (51) 
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gives 

5" = J dzq'il-qf. (52) 

This expression is only valid for q : ^ 1; given the derivation, it is a periodic 
function in q with period one. Also note that in anticipation of later results, I 
have replaced a factor of the length L by J dz. 

Physically, the computation of quantum corrections has lifted the degeneracy 
in q by generating a potential for q. In the full effective action, it helps to rescale 
the length in the z direction, introducing 

z' = \ll9Tz- (53) 
With this rescaling, the complete effective action at one-loop order is 

geff ^ + ^ Li S' , (54) 

9 

where 

S' = j dz' ((^)' + ^{l-lf^ ■ (55) 

The factor of follows on dimensional grounds, as at high temperature, T is 
the only natural mass scale in the problem. (The renormalization mass scale 
doesn't appear until next to leading order ^.) 

What is most interesting is that the which we had expected, becomes 

only a l/g. This is because our effective action only acquires a potential at one 
loop order. 

This begs the important question, why should this effective action be trusted 
to one loop order? What about the effects to higher loop order? The point is 
that in the new effective action, the relevant distance scale is not just l/T, but 
l/{gT): notice the factor of gT in the definition of the rescaled length z'. Thus 
for small g, any variations in the effective action occur over much larger distance 
scales than l/T. This is why our method of derivation — ignoring the variation 
of q{z) in the quantum action — - works. It does vary in space, but in weak 
coupling, this variation is very slow, and can be ignored. 

Having reduced the effective action to the above form, we merely want the 
"kink" which interpolates between the two vacua. While the general form of 
the kink is well known, in fact we only need the action. To compute the action, 
we note that the "energy" e for this system is conserved: 

(56) 

Here we view the spatial coordinate z' as a kind of time; saying the energy is 
conserved means that it is independent of z' . The energy is then a constant 
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of motion; with these boundary conditions, this energy vanishes. Zero energy 
impHes 

^^,(1-,)^ (57) 



Using this, 



S' = 2 J dz'q^{l-qf = 2^ dqq{\-q) = \. (58) 



3 



Putting everything together, 



S^ff ^ Ll^-, (59) 
3V6 <? 

which is the final result for two colors. 

This demonstrates that ci ^ 0; the interface tension vanishes classically, 
but is generated through quantum effects. The most interesting feature of the 
analysis is how the tunneling amplitude goes from the expected 1/g^ to just 1/ g. 
There are examples in string theory where tunneling amplitudes are not l/g^, 
but only \/g. These examples appear special to string theory, as the appearance 
of the coupling constant in this fashion is more or less natural. That is not the 
case here; the 1/g is really novel. The appearance of the inverse distance scale 
gT is reasonable, as the Debye screening mass in a thermal bath. 



3.2 Spins and Polyakov Loop Spins: Duality? 

Remember the behavior of a usual Ising magnet, in which spins Gi — ± interact 
through a coupling constant Jmag- The Hamiltonian is: 

Ti. = — Jmag '^i,n C i ■ <J i+h ■ (60) 

The sum is over all lattice sites, i, and nearest neighbors to «, fi. The spins align 
at low temperatures, and disorder at high temperatures. This is just because 
the partition function is Z ^ exp(— 7i/T). In magnets, the spin-spin coupling 
is more or less independent of temperature, so that with ~ Jmag/T in the 
exponential, ordering wins at low temperature, and loses at high temperature. 

The interface tension can easily be estimated. The simplest interface is to 
take all spins on the left hand side spin up, and all on the right hand side, spin 
down. If a is the lattice spacing, then the interface tension, defined as above, 
is J; by construction, its width is a. This very sharp interface is not the 
configuration of lowest energy, but the true interface tension is of order ~ J, 
with a width of order, a. Another example, more familiar to field theorists, is 
given by a scalar field with a double well potential |Q . 

Now consider an effective lagrangian for Polyakov loops. Over distances 
> l/T, the four-dimcnsional theory reduces to an effective theory of spins in 
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three spatial dimensions. For two colors, Polyakov loops are a type of Z{2) spin, 
with £ — ±. From the above S'^f^, 

'-^Polyakov ^ ■ (^'^) 

Now it is easy to understand why Polyakov loop spins order at high, instead 
of low, temperature. For magnets, the partition function involves exp{~ Jmag /T); 
for Polyakov loop spins, instead we have exp{—Jpoiyakov/T); but as Jpoiyakov ~ 
T^, in all the temperature dependence in the exponential is not l/T, as for or- 
dinary magnets, but T! 

This also leads to a natural conjecture of duality: that the temperature for 
Polyakov loop spins, and an ordinary magnet, are related as 

1 

Tpolyakov ^ 7^; • V^^J 
mag 

I have assumed that the variation of the gauge coupling constant with temper- 
ature can be neglected. 

This argument is extremely heuristic, and carries an important qualification. 
For ordinary spins, the lattice spacing is of course fixed. (This is true as well for a 
scalar field with double well potential [^.) From the derivation of the interface 
tension above, however, the width of the interface is the inverse Debye mass, 
^ l/{gT). Thus the argument fails in the limit of high temperature, since then 
the size of any single domain is becoming very large, ~ \/log(T)/T, as T — > oo. 
This doesn't contradict the conclusion of ordering at high temperature, since 
any single domain is, by definition, an ordered state. 

Now assume that the transition is of second order, as happens for two colors 
[ po| . Then the correlation length diverges at Tc] with Polyakov loop spins, it 
decreases as T increases from Tc- This divergence is determined as usual by 
scaling at a critical point. Even so, the underlying length scale which fixes the 
lattice spacing for effective Polyakov loop spins is fixed, set by a mass scale 
proportional to Tc, etc. Near Tc, we have implicitly made the assumption that 
the interface tension remains proportional to ~ T^. This cannot be true very 
near Tc, since for a second order transition, the interface tension must vanish 
at Tc. 

So is the argument of any use? Well, if d is the number of space-time 
dimensions, then simply on geometric grounds, Jpoiyakov ~ T'^"^. In three 
dimensions, then, Jpoiyakov ~ T; depending upon the value of Jpoiyakov/T, the 
system can still order above Tc, so there is no obvious contradiction. However, 
it does suggest that the width of the critical region is much narrower in four, as 
opposed to three, dimensions. For SU (2) gauge theories, this comparison is of 
interest in its own right. 

For a strongly first order transition, at first sight one might think that one 
should be able to directly check if Jpoiyakov ^ T"^^^. This is complicated by the 
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fact that near Tc, Z{N) states don't tunnel directly from one to another, but 
from one Z{N) state, to the symmetric vacuum, to another Z{N) state. 

3.3 Polyakov Loop Spins: Potts versus GKA 

The above analysis can be extended to more than two colors. For two colors, 
there is only one interface tension, between i = +1 and ^ = —1. For N colors, 
the vacuum is one of the A^th roots of unity, I = exp{27Tij /N), j = 1 . . . {N — 1). 
By charge conjugation, under which £ i* , the states j and N—j are equivalent. 
There are then about ~ N/2 distinct interface tensions. 

At any N, the smallest interface tension is between j — and j — 1. Defining 
j^e/z _ aiLj, then at next to leading order, 

ar = ^^^^4^ (1 - 12.9954...^ + ...), (63) 
3V37V 9{T) \ (4^)2 ; ' ^ ' 

where the running coupling constant g^{T) is defined using a modified MS 
scheme 

I remark that Boorstein and Kutasov [|4| argued that due to infrared di- 
vergences, (71 is not ~ 1/5, but one over the magnetic mass scale, cti ^ 1/ff^- 
While hardly conclusive, at least at next to leading order, there are no sign of 
infrared divergences. 

The interface tension from j — Q to j — k has been computed by Giovan- 
nangeli and Korthals Altes (GKA) The result is amazingly simple: 

k(N-k) 

"fc = N-1 ■ ^^^^ 

Now I construct an effective theory of discrete Z{N) spins. I forget about 
the factors of temperature which preoccupied me in the previous subsection; all 
that I am concerned with is the dependence on the distance between the Z{N) 
spins. This suggests what I term the GKA model. The spins at each site of the 
lattice are integers j, j = . . . {N — 1); the Hamiltonian is 

Ti-GKA = JgKA '^ji k {N ~ k) , k = \ji - ji+nlrnod N ■ (65) 

Tracing through the factors of A^, and holding g'^N fixed as ^ oo, the coupling 
constant Jgka > is of order one as ~> oo. 

The GKA model is in contrast to the Potts model, with Hamiltonian 

Ti-Potts = J ^ji 4o • (66) 

For the Potts model with J > 0, the energy is lowered if if two spins are equal, 
while if they differ — no matter by how much — the energy vanishes. 

For two and three colors, there is no difference between the GKA model and 
the Potts model. For example, consider the case of three colors. Then j = 1 is 
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equivalent to j = 2, so there is only one interface. That is, for the three roots 
of unity, any root is right next to the other two. 

The Potts model is known to be of first order for any number of states 
greater than, or equal to, three. For the GKA model, in mean field theory the 
transition is of first order for four colors ||3^; after all, the interaction between 
j = and j 2 is 4/3 that between j = and j = 1. Thus it would be very 
surprising if the GKA model wasn't also of first order when = 4. Further, 
Lattice simulations of SU (4) find a first order deconfining transition . 

As the number of colors increases, though, the Potts and GKA models be- 
come increasingly different. Whatever N is, in the Potts model any spin state 
interacts equally strongly with any other spin state. In the GKA model, at large 
N spins only interact significantly with those which are close in spin space. For 
example, ai ^ N, while aj ^ N"^ for j ^ N/2. It would be interesting to know 
the order of the phase transition in the GKA model at large N , both in mean 
field theory and numerically. 

This assumes that the interaction between Polyakov loop spins — computed 
in the limit of very high temperature — remains the same all of the way down 
to Tc- This certainly is wrong for a second order transition, but the question 
here is if it is first order. Thus: does the interface tension stay ^ j{N — j) 
(higher powers of ~ j'^{N — jY only make the more less like Potts), or become 
constant, independent of j? 



3.4 Z{N) Bubbles 

From the one-loop effective action, we can define a "potential" for q due to 
gluons, Vgi. I will be schematic, suppressing all inessential details, and taking 
two colors for now. From the computation of the one-loop effective action, it 
is clear that the assumption of Lt <C i was actually a matter of words; one 
obtains identically the same result in an infinite volume. Thus from S'' , I define 

V,i{q) - {q'{l~qf + f,) (67) 

/gT^ is proportional to the free energy of gluons at a temperature T . As noted 
above, this potential is only valid in the region < g < 1; it is periodic, with 
period one, outside of this region: 

Vgi{q+l) = Vgiiq) . (68) 

Obviously, q — Q and q — 1 are degenerate, 

V,m = V,i{l) . (69) 

This follows from the Z{2) symmetry of the pure glue theory. 
The quark contribution is computed similarly: 



V,k{q) - {2q^ - q^ + fqk) ; 



(70) 
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fqkT'^ is proportional to the quark contribution to the free energy. This expres- 
sion is only valid for > g > 1; else q is defined modulo one. Consequently, 

V,kiq + 2) = Vgkiq) . (71) 

This must be true for any potential, since q = and q = 2 both give i = +1. 
Moreover, while g = 1 is a an extremal point of the potential, it is a local 
maximum, not a minimum, with 

V.kil) > V.kiO) . (72) 

This shows how quarks violate the global Z{2) symmetry of the two color theory. 

To one loop order, the total potential for q is the sum of the gluon and 
quark contributions. With many (light) quark flavors, the total potential is like 
the quark contribution, with a maximum at q — 1. Dixit and Ogilvie |]38| first 
noticed that if the number of quark flavors isn't too large (or if the quarks are 
sufficiently heavy), q = I can be a local minimum; i.e., 9 = 1 is metastable. 

The above carries through for an arbitrary number of colors and flavors. If 
metastable states arise, they are necessarily a Z{N) state, with i a (nontrivial) 
A^th root of unity. They are termed "Z{N) bubbles" [|9[|38|. 



This all appears to be directly analogous to the usual problem of metastable 
vacua, but there is one important difference. For an ordinary potential, either 
in quantum mechanics or in field theory, one never worries about the zero of 
the potential, as that can be shifted at will. In the present case, however, the 
zero of the potential is physical, and gives the free energy of the stable vacuum. 
This is because the "potential" is multiplied by an overall factor of T*, and 
thcrmodynamically, derivatives with respect to the temperature matter. 

Thus there is no freedom to change the zero of the potential for q. For some 
Z{N) bubbles, if the potential at q 7^ is much higher than q = 0, it is well 
possible that the pressure in the bubble isn't positive, but negativel This was 
noticed first by Belyaev, Kogan, Semenoff, and Weiss js^. 
This is a complete disaster. I suggest a possible resolution. 
When we deal with g, we are in fact dealing with an angular variable. If we 
write the potential for q in terms of the thermal Wilson line, it is 

V,.rt{q) - {l~q? - r4tr((logL)2(l - (logL)^)) . (73) 

This form is correct in perturbation theory, where at each point in space, 
L(a;) is an element of SU (N). As will become clear in the next section, however, 
if we construct an effective theory for L, it no longer is an element of SU{N). 
Then this potential, for the purely angular part of L, is ill defined when its 
modulus vanishes. This ambiguity is easily cured by multiplying by an overall 
factor of the modulus: 

Vnon-pertm ^ (I^P + ■ ■ •) tr ((logL)2(l - (logL)^)) . (74) 
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This is rank conjecture: it is certainly a noK-perturbative modification of the 
potential. There is no reason to exclude terms of higher order in ~ 

This still does not solve the problem of the zero of the potential. I now 



assume further that the Polyakov Loop Model (PLM) applies |26, 



pO| , |31[ . Ignoring the angular variation in logL, the PLM potential is 

VpLMii) - T^b2\i\' + b.m^) . (75) 

The exact potential depends upon the number of colors and flavors, etc., but 
this is inessential here. For q — 0, the "usual" free energy is given by minimizing 
Vplm(^) with respect to £; with the above convention, the "mass" squared for £ 
is negative in the deconfined phase, 62 < 0, and positive in the confined phase, 
62 > 0. 

The complete potential is the sum of V„o„_pert(L) and Vplm(^)- At fixed 
£, as before any metastable points are elements of Z{N). The equation which 
determines t, however, is changed, as any metastable state has an action which 
acts like a positive mass term for t. Thus the expectation value of £ in a Z{N) 
bubble — even deep in the deconfined phase — has £ < 1, not £ = 1. Before a 
Z{N) bubble develops negative pressure, £ = 0, with zero pressure in the PLM. 
More likely, Z(N) bubbles become unstable in the £ direction before £ — Q. 

This could be tested on the Lattice. Compute in a theory in which the 
splitting between the true vacuum and the metastable state is small; (dynamical) 
heavy quarks will do. Then the expectation value of I should be smaller in the 
metastable state than in the stable vacuum. This holds regardless of the question 
of renormalization discussed in the next section. 

If true, all perturbative calculations of the lifetime of a metastable Z{N) 
bubble are wrong [Q. At best, they are a upper bound on the true lifetime, 
which is somewhat useless. On the other hand, it resurrects the possibility that 
Z{N) bubbles — which spontaneously violate CP symmetry — might have 
appeared in the early Universe. 



4 Renormalization of the Wilson Line 

The sections following this deal with mean field theory for the thermal Wilson 
line. Implicitly, this assumes that it is possible to go from the bare Wilson line, 
as measured on the Lattice, to the renormalized quantity. How to do this on 
the Lattice is presently an unsolved problem; in this section I review what is 
known |§. 

In a pure gauge theory, the expectation value of a closed Wilson loop, of 
length L and area A, is 

(tr'P explig / = exp (-mo-L - ctA) . (76) 
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The string tension, a, is nonzero in the confined phase, T < Tc, and vanishes in 
the deconfined phase, T > Tc- 

The concern here is not with the term proportional to the area of the loop, 
but with the length. This is a type of mass renormalization for an infinitely 
heavy quark. For example, it is easy to compute this to lowest order in per- 
turbation theory. We are interested in an ultraviolet divergent term, so over 
short distances, it suffices to assume that the loop is straight. For the sake of 
discussion, assume that the loop runs in the time direction. (New divergences 
arise when there are cusps in the loop; these divergences can also be computed 
perturbatively, by a similar analysis j40|.) Then to lowest order, there is a 
contribution 

~ -g^ ( dn f dT2 Ao{x,T,)Ao{x,T2)) ^ f d^k^. (77) 

Jo Jo ^ J 

The integral is nominally divergent, but with either dimensional or Pauli-Villars 
regularization, the divergence vanishes l ^of . This cancellation is somewhat triv- 
ial at one loop order, arising from having three powers of momentum upstairs, 
and two powers downstairs. 

Thus one would expect divergences to arise at ~ g^, which are found. How- 
ever, for closed Wilson loops, all such divergences can be absorbed into charge 
renormalization |4^. This is a notable result: in a quantum field theory, gen- 
erally the renormalization of any composite operator requires the calculation of 
its mixing with all other operators with the same mass dimension and symme- 
tries. Like the action itself, however, the Wilson loop has a privileged status; it 
doesn't mix with any other operator. 

As noted first by Polyakov |40j, this result can be understood on the basis 
of reparametrization invariance for the Wilson loop. We parametrize the loop 
as a curve x^{s)^ where s is the length along the path. Then with x — dx'^/ds, 
the term 

J Vi^ ds (78) 

is invariant under s s'{s). Generally, physics shouldn't depend upon how we 
label path length along the curve. Because it has dimensions of length, however, 
the coefficient of this term must have dimensions of mass. With dimensional 
regularization, there is no such mass scale. (The renormalization mass scale 
only enters to ensure the proper running of the coupling constant.) A term 
which has no mass dimension is 

x = d^x^ /ds^. This is not reparametrization invariant, though, and so does 
arise with dimensional regularization. 
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On the other hand, assume that the regularization scheme does introduce a 
mass scale. On the Lattice, this is the inverse lattice spacing, ^ 1/a. Then a 
term proportional to the length, L, does appear pO| , 

- -g' I- (80) 

At nonzero temperature, L/a = Nt, the number of lattice steps in the time 
direction. Clearly, this is the first term in an infinite series in the coupling 
constant, g. 

How to deal with the power divergences generated by the Lattice is at present 
an unsolved problem. Since in the continuum there are neither logarithmic nor 
even finite terms to worry about, this appears to be a technical, albeit important, 
problem to solve. 

Why is this important? In the confined phase, this constant is not of any 
particular consequence, as the Wilson loop is dominated by the string tension. 
In the high temperature phase, however, the trace of the thermal Wilson line is 
the order parameter for the phase transition. It would be peculiar if a precise 
physical definition did not exist. Any composite operator requires a condition to 
fix its renormalized value; for the thermal Wilson line, the natural prescription 
is that Polyakov loops approach one as T ^ oo. 

To compute the leading perturbative correction to the thermal Wilson line, 
it is necessary to include effects from the Debye mass, rriD ^ gT. Replacing the 
bare propagator for Aq, by 1/(P + 

T J +m1 

This divergent integral can be computed using either dimensional or Pauli- 
Villars regularization. Or, one can just subtract the integral with m^) = 0: 

This was first demonstrated by Cava and Jengo That is, while the leading 
term is negative in the bare theory, it is positive after regularization. This 
change in sign is unremarkable, as the sign of a renormalized operator is not 
preserved under regularization. 

This appears to indicate that the renormalized thermal Wilson line is not a 
unimodular matrix. One concern is that any quantity g^ really arises from a 
two-loop graph, ^ g^, times an infrared singular piece ~ l/rriD ~ Thus it 
is necessary to ensure that the above is the only infrared singular term at this 
order. 

A different calculation was performed by Korthals Altes He computed 
the one loop corrections to the thermal Wilson line in a background Aq field. 



(81) 



Notes on DcconGncmcnt 



22 



The method is identical to that used in sec. II to compute the interface tension. 
ClassicaUy, the thermal Wilson line is a special unitary matrix. Korthals Altes 
finds that the one-loop corrections to the thermal Wilson line are not only 
infinite (!), but generate a matrix which is neither unitary nor special. On the 
other hand, all Polyakov loops are finite. 

Thus even in the continuum, the renormalization of the thermal Wilson line, 
and Polyakov loops, remains an unsolved problem. 

A way of measuring renormalized Polyakov loops on the Lattice has been 
proposed by Zantow et al. They compute only two-point functions of 

the Polyakov loop. At short distances, the static potential can be computed 
perturbatively, which allows one to extract the renormalized Polyakov loop. 



5 Deconfining Transition for Two, Three, and 
Four Colors 

5.1 Polyakov Loops and Quarkless Baryons 

So far I have been concerned with the (pure glue) theory at very high tempera- 
tures. Now I turn to the question of its behavior near the critical temperature. 
I review Lattice results on the order of the phase transition for two ||2^, three 
[ pT| , and four ||l7|| colors, and then ask what constraints it places on the mean 
field theory for Polyakov loops. 

Up to this point, I have only considered the trace of the thermal Wilson 
line in the fundamental represenation, which is the the Polyakov loop t. By 
a local gauge transformation, at each point in space one can diagonalize the 
thermal Wilson line. These eigenvalues are gauge invariant, so since L(x) is an 
SU (N) matrix, at each point there are N — 1 independent degrees of freedom. 
Another — 1 degrees of freedom are given by the trace of powers of L, trL-' , 
j = l...iN-l). 

Under a global Z{N) transformation, the "usual" Polyakov loop transforms 
as a field with charge-one, eq. (^; thus I relabel it £i. Polyakov loops with 
higher Z{N) charge are easy to construct. I introduce the traceless part of L: 

L = L-£il. (83) 

Then I define the charge-two Polyakov loop to be 

^2 = ^trL2 = ^trL^ - ^ (trL)^ , (84) 

where 

l2 -> e2'*£2 , (85) 
with (j) as in eq. (|^). There are two operators with charge-two, £2 and £\. 
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For example, consider two colors, and the parametrization of the thermal 
Wilson line in the strict perturbative regime, (24). The charge-one loop is 
£i = cos(7rg), eq. (|2^), while the charge- two loop is £2 = — sin^(7r(7). At high 



temperature, where g = 0, 1, {£1) ±1, while {£2) 0. 

I note that Polyakov loops of charge-two and beyond are related to the trace 
of the thermal Wilson line in higher SU{N) representations. For two colors, in 
perturbation theory the trace of the Wilson line in the adjoint representation is 
trCLadj) = 1 + 2^2 • So far, though, I haven't found this particularly useful. 

Continuing on, 

4 = ^trL3 (86) 

has charge three under the global Z{N) symmetry. Other charge-three operators 
are £1 and £i£2- The construction of operators with higher Z{N) charge proceeds 
similarly. For example, operators with charge four, independent of the singlet 
part, are given by trL^ and (trL^)^, etc. 

I stress that both the expectation values, and correlations functions of, 
Polyakov loops of arbitrary charge are well worth measuring on the Lattice. 
When the Z{N) symmetry is spontaneously broken at T > Tc, Polyakov loops 
£j with charge j = 1 . . . {N — 1) all acquire nonzero expectation values. As 
Polyakov loops of charge-two and beyond are constructed from the traceless 
part of the thermal Wilson line, they aren't that interesting at high tempera- 
ture; as T —^ 00, their expectation values are proportional to nonzero powers of 
g^, times powers of T to make up the mass dimension. Near Tc, however, there 
is nothing general which can be said about their expectation values. 

I will assume that for a second order deconfining transition, the only critical 
field is the charge-one loop Even so, Polyakov loops of charge- two and 
beyond will certainly affect non-universal behavior. One notable example is 
the Polyakov Loops Model |2^, H^, 30, which conjectures a 



relationship between the expectation value of Polyakov loops and the pressure. 
The original model assumed that only the charge-one loop mattered, but I no 
longer see why the expectation values of higher-charge loops are not important 
as well. 

There are certain Polyakov loops which have a privileged status: these are 
those with charge- TV. As they are neutral under Z{N), their expectation values 
are nonzero at any temperature. I term such operators quarkless baryons. 

In QCD, a baryon is N quarks tied together through an antisymmetric tensor 
in color space. One can also consider a more general object, a baryon "junction" 
This is an antisymmetric color tensor, with N Wilson lines coming out of 
it. Putting quarks at the end of each line gives the usual QCD baryon, since 
in the confined phase, Wilson lines are short, on the order of ^ l/y^, where a 
is the string tension. While directly related to QCD baryons, without quarks, 
baryon junctions are not gauge invariant: only junction anti-junction pairs are. 

In contrast, all quarkless baryons are gauge invariant. The simplest quarkless 
baryon is £^ . In mean field theory, this is zero in the confined phase, and nonzero 
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above. In the full quantum theory, {i^) 7^ at all T. While there is not good 
Lattice data on this expectation value, I assume that it is small below Tc, and 
large above, but this is just a guess. Since junction anti-junction pairs involve 
N Wilson lines, they are directly related to the operators for quarkless baryons, 
etc.. 

5.2 Effective Theories for Polyakov Loops 

I next turn to the construction of an effective theory for the thermal Wilson line. 
Remember how this proceeds with an Ising model on a lattice. While the value 
of the spin on each site is ±1, after an effective spin is computed by averaging 
over a domain of fixed size, the result effective spin is a continuous variable, 
4>{x). The effective theory is just the usual t/)^ theory. By the renormalization 
group, it is in the same universality class as the original Ising model. 

The analogous proceedure can be carried through for the thermal Wilson 
line. The effective thermal Wilson line, constructed by a gauge invariant 
average over a domain of some size, is not an SU{N) matrix, but has more 
degrees of freedom. I then consider all Polyakov loops, from charge-one up 
to charge-iV. Of course there is no reason to stop there, but presumably the 
number of effective fields which really matters is limited. 

In an effective Lagrangian, the first thing to ask about are the mass terms: 

C'ff = m2|£i|2 + to2|^2|2 + ... (87) 

The simplest assumption is that for T > Tc, condensation is driven by a negative 
mass term for the charge-one loop: 

TO? < , T > Tc , TO? > , T < Tc, (88) 

and that the masses for all higher loops are positive at all temperatures, 

ml>0 , TO^ > . . . (89) 

If so, then the charge-one loop controls the critical behavior 

There is good reason why one expects that condensation is driven by that 
of the charge-one loop, and not by loops with higher charge. If the mass for the 
charge-one loop is negative, the favored vacuum is given by maximizing |trLp. 
After a global gauge rotation, we can always choose the expectation value of L 
to be a diagonal matrix. If L were a U{N), instead of an SU{N) matrix, then 
|trLp is maximized when L is a constant phase times the unit matrix. This 
remains true when L is a SU{N) matrix; for it to be a unit matrix, however, it 
must be an element of the center of the gauge group, 

(L> = ^0 exp 1 . (90) 



with a Z{N) phase, = ^irj/N, j ^l...{N -I). 
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For this particular expectation value, the vacuum does not spontaneously 
break the (global) SU{N) symmetry above Tc. This accords with naive expec- 
tation: the high temperature vacuum is not in a Higgs phase. The possibility 
of having an expectation value which doesn't break SU{N) is special to a field 
in the adjoint, as opposed to the fundamental, representation. 

On the other hand, assume that mf > 0, and m| < in the deconfined 
phase, so that symmetry breaking is driven by condensation of the charge- 
two, instead of the charge-one, loop. Then the vacuum is given by maximizing 
jtrL^P; this means that the expectation value of is an element of the center. 
But if so, besides the SU{N) invariant vacuum, there are also vacua which are 
only invariant under SU{N — 1). At present, there is no evidence to suggest 
that the high temperature vacuum is one where SU{N) spontaneously breaks 
to SU{N-1). 

The above description can be extended beyond mean field theory, at least if 
the deconfining transition is of second order. A transition driven by the charge- 
one loop is one where rnf —> at Tc; one driven by the charge- two loop is where 
m| — > at Tc, etc. To be precise, for a transition driven by a charge-fc loop, 
both its mass, and that of the charge-(A^ — k) loop, vanish at Tc. 

The masses for Polyakov loops of all charges can be directly measured on 
the Lattice. Even for the charge-one loop, data near Tc is, at present, limited 
[|[ |22| . There is also some data for higher charge loops for three colors in 2 -I- 1 
dimensions [ ^ . 

Given this (crucial!) assumption about the masses, I next turn to the order 
of the phase transition for a small number of colors. 

5.3 Two Colors: Second Order, and Quarkless Baryons 

For two colors, all Polyakov loops are real. For the charge-one loop, I take the 
potential 



with a positive quartic coupling, Ai > 0. Of course higher powers in £i are 
also possible. Near a second order phase transition, however, the most relevant 
operators, with the fewest powers of ii, dominate. 

Invariant terms involving the charge-two loop include 



plus terms ~ ^ ^^c. All powers of £2 are allowed because it is a singlet. 
The potential which mixes the charge-one and charge-two loops starts out as 




(91) 



V2 ^ h£2 + - 




(92) 



V, 



^£l£2 + -- - ■ 



(93) 



plus many other terms; this has the lowest mass dimension. 
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There is extensive Lattice data on the nature of the deconfining phase transi- 



tion 1 20 . Especially from the work of Engels et al, it appears that the transition 
is of second order. To wit, the critical exponents are within ^ 1% of the values 
expected for the Ising model |^ . 

There is a surprise, however. As first stressed by Damgaard |20[ |, the expec- 
tation value of the Polyakov loop in the adjoint representation is an approximate 
order parameter. In perturbation theory, the adjoint Polyakov loop is 1 -I- 2£2, 
so from the Lattice data, the expectation value of £2 presumably jumps at T^. 

From the terms above, the expectation value of the charge-two loop is 

(4) = _ '1±^ . (94) 

To explain the jump in (£2) about Tc, there are then two possibilities. If the 
charge-two loop is heavy, then the coupling constant of the charge-two loop to 
the charge-one loop, ^, must be large. This means that changes in the density 
of the charge-two loop is driven by condensation of the charge-one loop. 

The other possibility is that h and ^ are not especially large, but that the 
charge-two loop becomes light near Tc- The latter doesn't violate universality, 
as long as the charge-two loop isn't massless at Tc- 

There is no lattice data on (il). I presume that as suggested by mean field 
theory, quarkless baryons are rare below Tc, and common above. This is 
seperate from the changes in £2- 



5.4 Three Colors 

5.4.1 A "Nearly" Second Order Transition 

In an asymptotically free gauge theory, it is natural to form the ratio of the 
true pressure to that of an ideal gas. In principle, positivity of entropy does not 
require this ratio to be less than one. In practice, Lattice data with improved 
actions finds that this ratio is less than one at all temperatures . 

Numerical simulations find that for three colors, the deconfining transition 
is of first order , in agreement with general arguments ||] . For a first order 
transition, the pressure is continuous at Tc, but the energy density jumps. Thus 
consider the ratio of the jump in the energy density to that of an ideal gas. 

This ratio is not bounded by one. To illustrate this, consider a bag model. 
Above Tc, the pressure is that of an ideal gas, minus a bag constant, b: 

Pbag = CqT'^ - b . (95) 

The pressure is assumed to vanish below Tc, with Tc fixed by pbag — 0. This 
bag model does not describe the Lattice data near Tc, but is a useful construct. 
In the bag model, the ratio of energies is: 
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In contrast, Lattice data appears to give a result which is much smaller: 

Se . 1 , , 

{Lattice ^ - \y ' } 

The mass of the charge-one loop has also been measured on the Lattice |^ . 
It goes from rrii/T ~ 2.5 at ^ 2Tc, down to rrii/T ~ .25 at ^ r+. This decrease 
in the screening mass, apparently by a factor of ten, strongly suggests the in fact 
the transition is even weaker than the above comparison with the bag model 
suggests. Instead of weakly first order, I prefer to call the deconfining transition 
nearly second order. 

(It is necessary to measure correlation functions of Polyakov loops to see 
this decrease. Masses measured from other operators, such as plaquettes, do 
not decrease dramatically about Tc This implies that the mixing between 
Polyakov loops and and plaquette operators are small. This small mixing is 
found in related problems Q.) 

An effective theory cannot explain why the deconfining transition is weakly 
first order; it merely requires that certain coupling constants are small. 



5.4.2 Polyakov Loops with Charge One and Minus One 

For three colors, I consider Polyakov loops with charge one, two and three. 

There is no data on the expectation values for the quarkless baryons, £f, 
£i£2, and £3. I presume that as indicated by mean field theory, {£1) is small 
below Tc, and large above. It would be interesting to knowhow the density of 
the other quarkless baryons, £1^2 and £3, change about Tc. Regardless of the 
renormalization issues discussed in sec. (4), changes in these expectation values 
are presumably physical. 

Thus I concentrate on the interaction between the charge-one and the charge- 
two loops. Remember that for three (or more) colors, the ij^s are all complex 
valued fields. The potential for charge-one loops is dictated by the global Z{3) 
symmetry: 

Vi = ml + {el + (eif) + Ai (|£i|2)' . (98) 

The notable feature is the appearance of a cubic term, which necessarily ensures 
a first order transition 

The charge-two loop has charge minus one under Z{3), so its potential is the 
same, albeit with different masses and coupling constants: 

V2 - mj + «2 {el + (^2)') + A2 {m'f ■ (99) 

There are many terms by which the charge-one and charge-two fields can 
mix. The most important is that with the smallest mass dimension: 



(100) 
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In terms of the original thermal Wilson line, this term is ~ (trL)(trL^), etc. 

If the charge- two field is heavy near Tc, we can integrate it out. While it 
may be a mess to do so analytically, any resulting potential, involving only li, 
must still respect the overall Z{3) symmetry. This produces a potential identical 
in form to Vi, but with different values for the mass and coupling constants. 
A weakly first order requires that the effective cubic coupling in the resulting 
eflfective theory is small, Ki <C 1. 

If the charge- two loop becomes light near T^, and if it mixes strongly with 
the charge-one loop through a large coupling constant ^, then its effects cannot 
be neglected. 

There is another possibility. The mass and quartic terms in the potentials 

are invariant not just under Z{3), but imder a global U{1) symmetry. Assume 
that the charge-two field is always heavy. Then all terms in both potentials 
are invariant under a global U{1), with the exception of the cubic terms, with 
couplings Ki and K2, and the mixing term between £1 and £2, with coupling ^. 
Thus perhaps all terms invariant under Z{3), but not U(l), are small. That is, 
the heavy charge-two loop mixes weakly with the charge-one loop. This doesn't 
explain why all Z{3) couplings are small, but hints at a more general principle. 

It will be interesting to see what detailed numerical studies on the Lattice 
tell us. 



5.5 Four Colors: First Order from Charge-Two Loops 

For four colors, I consider just the charge-one and charge-two loops. Under 
Z{4:), £1 i£i, so the potential for the charge-one field alone is 

2 

Vi = ^ Kip + Ai {\£,ff + K, (if + {£1)^) . (101) 

The term ^ Ai is 0{2) invariant, while that ^ Ki is only invariant under Z{4). 

Under Z{4:), £2 —£2, so the potential for the charge-two field by itself is 
just like that for the charge-one field: 

V2 = +ml \£2f + A2 {\e2ff + K2 {£t + {£*2)^) . (102) 
The allowed terms which mix the two fields are: 

Vrmx = Cl (^;^? + h{£lf) + C2 {e2£l+£*2iilf) ■ (103) 

Unlike two colors, a term linear in £2 is not allowed by the Z(4) symmetry. 
For (2=0, and assuming that the charge- two field reamins heavy about T^, I2 
can be integrated out to give: 

~ -^(Kil^)^ (104) 

7712 
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One can convince oneself that this term is necessarily negative. That is, the 
quartic coupling for the charge-one field is shifted downward: 

Ai EE Ai - . (105) 
mi 

The same holds if both and ^2 are nonzero: integrating out the charge-two 
field generates corrections to the quartic coupling constants of the charge-one 
field which are uniformly negative, shifting them downwards. 

The transition for four colors appears to be of first order ]l7t . One explana- 
tion for this is that the original coupling constants Ai and A2 are positive, but 
after integrating out the charge-two field, they become negative, and thus drive 
the transition first order. 

The analysis for higher numbers of colors is then immediate. The leading 
term which couples a charge-j loop to the charge-one loop is 

V™„ = ^(^£*£\ + ijiiiy) . (106) 

Assuming that the charge-j loop is heavy about Tc, we can integrate it out, 
which produces a term in the potential for the charge-one loop ~ (j^ip)-'. For 
charges greater than two, this is less relevant (has smaller mass dimension) than 
the quartic terms expected to dominate. 

Thus if the deconfining transition is of first order for more than four colors, 
in the present language it is uniquely due to how the coupling between charge- 
two and charge-one loops affects the effective quartic coupling constant for the 
charge-one loop. Polyakov loops with charge greater than two do not affect the 
order of the transition. 

6 A Parting Comment 

While it is true that, in equilibrium, all thermodynamic quantities follow from 
the pressure, experience with the perturbative calculation of processes near 
equilibrium — such as transport coefficients, real photon production, etc. — 
teaches us they often depend on the details of equilibrium correlation functions. 
Thus regardless of theoretical prejudice, such as the Polyakov Loops Model, it is 
important to measure as many gauge invariant correlation functions as possible. 
There is an astounding amount of superb data which is pouring out of RHIC 
Many features, including the change in the spectrum of "hard" particles, 
details of Hanbury-Brown-Twiss interferometry, chemical composition, etc., ap- 
pear to defy explanation by any conventional mechanisms js^. We may very 
well need a detailed understanding of the theory, near Tc and above, in order 
to sort out these amazing results. 
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